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I . INTRODUCTION 

The algorithm proposed in this report for smoothing 
3-D data of NUWES uses sequential differences to screen the 
data for "wild" data values (outliers) and a 7-point least- 
squares procedure to perform the data smoothing. 

It is assumed that the actual relationship between 
each coordinate (x,y,z) of the path of the vehicle being tracked 
and time can be expressed as a low order polynomial P^(t) 
where k is the order (degree) of the polynomial. An observed 
value can then, for example, be expressed as 



x . 

1 



= P, (t. ) 
k 1 



+ n . 
1 



where n^ is the noise component resulting from inaccuracy in 
measuring the x-component of the position of the vehicle at 
time t^. Examination of 3-D data from a torpedo run at NUWES, 
supported by consideration of the relationship between 
vehicular maneuvering capabilities and data rate, suggests that 
polynomials of no higher order than K = 3 (cubics) need be 
considered for 7 consecutive point data segments. The selec- 
tion of the appropriate order of polynomial and the smoothing 
are incorporated in STEP 4 of the algorithm and discussed in 
Section III D. 

Performance of least-squares smoothing can be stand- 
ardized and computational requirements considerably reduced when 
7-point data segments contain no missing values. This neces- 
sitates provision of temporary values to fill in gaps in the 
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data. Such temporary values are also useful in the screening 
process to identify outliers. Provision of temporary values 
for missing data values is treated in STEP 1 of the algorithm 
and discussed in Section III A. 

When a component of a vehicular path is a polynomial 
of order k _< 3 , then the fourth order sequential differences 
of the observational data will consist exclusively of the noise 
components in the data, (Ref. 2) . Since outliers can be at- 
tributed to noise it would appear reasonable to use fourth 
order sequential differences to screen for them. This is in- 
corporated in STEP 2 and discussed in Section III B. This 
step calls for replacement of outliers by temporary values as in 
the treatment of missing values in STEP 1. 

The temporary values to fill in missing points and 
replace outliers were selected to be consistent with the closest 
observed data value before and the closest after the point in 
question. This could, and actually does in some cases, produce 
temporary values that are inconsistent with the polynomial that 
best fits the observed values in the data segment. It is pro- 
posed in STEP 3 that smoothing be performed as in STEP 4 for 
7-point data segments centered at the temporary values and that 
these smoothed values be treated as observed values in smoothing 
at observed data values. This is a direct attempt to minimize 
the effects of temporary values (.either missing points or 
outliers) on the results of data smoothing. This is discussed 
in Section III C. 
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The presence of two or more missing points and/or out- 
liers (i.e., temporary values) within the same 7-point data 
segment, or, more specifically, within three data points of 
each other, requires special treatment. Some initial efforts 
in this direction are included in STEP 5 and discussed in 
Section III E. The general motivation for this step is the 
derivation of the maximum information on a vehicular path that 
can be obtained from the data even though that information is 
degraded by missing points and outliers. The possibility of 
providing some measure of the "quality" of the smoothing at 
each data point for the guidance of potential users is discussed 
in Section III F. 
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II. A DATA SMOOTHING ALGORITH 



The following Algorithm using fourth order sequential 
differences and least-squares smoothing of 7-point data segments 
is proposed for smoothing of 3-D data at NUWES. 



STEP 1 


Preliminary treatment of missing data points. 


(a) 


A missing data value at time t^ where the adjacent 
observed values x^_^ an< ^ x i+i are ava il a ble is 
assigned the temporary value 

x ? = (x i-i + x 1+ i )/2 • 


(b) 


Adjacent missing values at times t^ and t^ + ^ , where 
the values and x i+2 are °^ serve< ^» are assigned 

the temporary values 

x* = (2x. i + x. ,»)/3 and 

l l-l i+2 

x l+l * (x i-l + 2x i+2 !/3 ‘ 


(c) 


Occurrence of more than two consecutive missing values 
are designated for Special Treatment in Step 5. 


STEP 2 


Identify and temporarily replace Outlier data points. 


(a) 


Calculate fourth order sequential differences d^.'s • 
(Ref. 2) 


(b) 


Any value d^ which exceeds a specified level d in 
magnitude is a potential outlier. The signature in 
fourth order differences of an outlier is a large value 
d^ with somewhat smaller values and ^4i+l 
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